The London Study of Vortex States in a Superconducting Film Due to a Magnetic Dot 
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Here we report a study of vortex states in a thin superconducting film with a magnetic dot grown 
upon it by means of a method based on London-Maxwell equations. Vortices with single quantum 
flux ($o = hc/2e), giant vortices ( vortices with multiple flux ) as well as antivortices ( vortices 
with negative vorticity ) are taken into consideration. It turns out that giant vortices occur, when 
the dot's size is sufficiently smaller ( R < 4.5£) than the effective penetration depth A. In the case 
of a dot with sufficiently large size ( R > 6£ ), the vortices with single quantum flux dominate the 
vortex states. Their geometrical patterns are predicted up to seven vortices. Our calculations do 
not show the spontaneous appearance of antivortices. 



PACS Number(s): 74.25.Dw, 74.25.Ha, 74.25. Qt, 74.78, 



Type II superconductors are used in a wide variety of 
technological applications due to their high critical cur- 
rents and fields [1]. In these superconductors, vortices 
appear when the magnetic field exceeds the first critical 
field H c \. Under external current or field, vortices move, 
which causes the superconductor to switch to a resistive 
state. As a result, the system loses its superconductivity. 
Because of this, vortex pinning is quite important in ap- 
plications of type II superconductors. One of the ways to 
pin vortices is to use the magnetic subsytems with either 
out-of-plane or in-plane magnetization. These subsys- 
tems are capable of trapping both vortices and antivor- 
tices, depending on the orientation and strength of their 
magnetization. The aforementioned systems are not only 
important for technological applications, such as devices 
that can be tuned by weak magnetic fields, but also offer 
rich physical effects which are not observed in the individ- 
ual parts. Some of these effects were predicted elsewhere 
[2-7]. 

In the recent decade, magnetic dots growing on top 
of SC films have been extensively studied both exper- 
imentally [8] and theoretically [9,10]. In experimental 
studies, magnetic dots with in-plane magnetization are 
fabricated from Co, Ni, Fe, Gd-Co and Sm-Co alloys, 
whereas, for the dots with magnetization perpendicular 
to the plane, Co/Pt multilayers are used [11] . These stud- 
ies report commensurability effects on transport proper- 
ties, which comfirm that the dots create and pin vortices. 
On theoretical side, several realizations of the aforemen- 
tioned systems are analyzed through the Landau-Gizburg 
framework and London theory. In these works, the au- 
thors investigated the conditions for vortices to appear 
and calculated their geometric configurations in equilib- 
rium. 

Recently, Priour et al. studied the vortex states of a SC 



film with a magnetic dot array grown upon on it, through 
Ginzburg-Landau theory, and found several different con- 
figurations of vortices in the SC film [12]. Similar study 
was also done by Peeters et al., but in the presence of a 
single dot [13,14]. One of the interesting results that is 
reported by these works is, when the dot's size is on the 
order of a few coherence lengths £, vortices with multi- 
ple flux quanta ( the giant vortices ) appear. The giant 
vortices are previously studied in the Landau-Ginzburg 
framework [15]. The stable antivortex (AV) ( vortices 
with vorticity opposite to vorticity of those trapped un- 
der the dot ) states were reported in similar systems [16]. 
The most interestingly, Kanda et. al. recently reported 
the experimental evidence of the GV states and the other 
rearrangements of vortices in a mesoscopic SC disc, using 
the multiple-tunnel-junction method [17]. 

Earlier, we studied the vortex states due to a ferromag- 
netic (FM) dot with out-of-plane ( perpendicular to the 
film ) magnetization on top of a SC thin film [19]. In that 
work, we limited ourselves only to the vortices trapped 
under the dot, and calculated the equilibrium configura- 
tions, up to three vortices. However, in a more realis- 
tic picture, both vortices under the dot and antivortices 
outside the dot's boundaries can appear spontaneously, 
because total flux due to the magnetic dot over the en- 
tire infinite SC film is zero. Actually, this problem was 
studied for magnetic dipoles with out-of-plane magneti- 
zation [20]. Another interesting case is giant vortex (GV) 
states with multiple unit flux $o. In usual circumstances, 
it is hard to get vortices with multiple flux quanta, be- 
cause their energy grows as the square of their vorticity. 
However, this energy cost can be overcome by the dots 
with radii of sufficiently less than effective penetration 
depth A = X 2 /d sc [18], where A is the London penetra- 
tion depth and d sc is the thickness of the SC film, and 
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with sufficiently large magnetization. 

In this article, we focus on the complete picture to- 
gether with antivortices and giant vortices, and aim to 
get analytical insight for the spontaneous vortex, antivor- 
tcx and giant vortex states. For this purpose, we pursue 
the method based on London-Maxwell equations, which 
is fully developed elsewhere [19]. Though the London ap- 
proach works well at the high n = A/£ limit, it perfectly 
serves our purpose and enables us to do analytical calcu- 
lations. In our analytical and numerical calculations, we 
show that there are three regimes for the vortex states 
according to the dot's size. Namely, when the dot's size 
is sufficently small, the GV states are dominant. If the 
dot's size is sufficently large, we see only vortex states 
with single quantum flux (SQF). There is also an inter- 
mediate region in which both GV and SQF states appear. 
We calculate the geometrical configurations of SQF states 
in equilibrium up to seven vortices. In addition, our cal- 
culations do not show any AV states. At this point, we 
need to remark that the London framework might seem 
unreliable to some readers, when the ratio of the dot's size 
to the coherence length, £, is in the order of unity. How- 
ever, one might expect that the vortices in the London 
framework, usually prefer staying in SQF states rather 
than GV states, since the London theory treats the vor- 
tex core as a point. Our approach then underestimates 
the dot's size, when the GV states appear. When the 
vortex cores are treated more accurately, as in Landau- 
Ginzburg framework, one might expect that the dot's 
size would be rather larger than that we estimate in this 
work. Nevertheless, the physics and qualitative result is 
still captured correctly in the London theory, although 
the results might be quantitatively different. 

In the following section, we briefly introduce the 
method we follow in this article. Next, we study the 
giant vortex states with and without antivortices. The 
third section is devoted to a discussion on the SQF vor- 
tex states with and without AV states. The last section 
consists of the discussions and conclusions. 



I. METHOD 



In order to attack this problem, we consider a thin cir- 
cular magnetic dot of radius R with magnetization per- 
pendicular to its plane at z = d, placed upon a thin SC 
him which both lay on the x-y plane z = (sec Fig. 1). 
The dot's magnetization reads 



m = mQ(R — r)S(z — d)z. 



(1) 
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SUPERCONDUCTING FILM 

FIG. 1. Magnetic dot on a superconducting film. 

The energy of the system in the presence of many vor- 
tices with arbitrary vorticities rn is described elsewhere 
[19]. The appearance of N vortices with arbitrary posi- 
tions Yi in the system changes the energy by an amount: 
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Here e v = £o m (A/£) is the self energy of a vortex with- 
out a magnetic dot, e = <!>§/(167r 2 A), e vv is the vortex- 
vortex interaction, and e mv is the vortex-magnetic dot 
interaction. The vortcx-vortcx interaction reads 
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where 



|rj — Ti\ , Hq(x) and Y (x) are the Struvc 



function of the zeroth order and the modified Bessel func- 
tion of the second kind of the zeroth order, respectively 
[21]. The asymptotics of Eq.(3) are 
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Note that in our further calculations, core energy e core <~ 
0.809eo is renormalized to A in logarithmic terms. The 
vortex-magnetization interaction energy e mv is repre- 
sented by 
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When d ^ and is significantly smaller than i?, the cor- 
rection to the e vm is on the order of d/R. Throughout 
the paper, we stick to d « R approximation to limit 
the number of parameters in our analytical calculations, 
which does not change the qualitative results. In our 
numerical calculations, however, we take d ^ to make 
sure that numerical integrations converge. The asymp- 
totics at d = are 
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Additionally, under the assumption that the magnetiza- 
tion does not change direction due to its interaction with 
vortices, the self-magnetization energy is irrelevant to our 
calculations. Because of this, we do not take this energy 
into account. In the rest of the paper, An denotes the 
energy of N vortices with single quantum flux (L = N 
state), while represents the energy of a giant vortex 
with vorticity N ( L g = N state ). A 9 ^ is used for a 
GV with N vorticity and N SQF antivortex around it, 
whereas A N is used for TV SQF states along with N AV 
states (L a = N state). 



II. THE GIANT VORTICES 

The self energy of a GV with vorticity n is proportional 
to its vorticity's square n 2 , whereas the self energies of n 
vortices with SQF is just proportional to n. Because of 
this, GVs are unstable in usual superconductors. Even if 
they pop up, they decay into vortex states with vorticity 
n=l. To get an idea about when GV states might occur, 
let us consider the following scenario in which SQF vortex 
states are trapped in a small region through a sufficiently 
small-sized magnetic dot so that they cannot go to fur- 
ther distances. In order for SQF states to be stable, the 
system requires more gain in energy, since energy grows 
due to the repulsion between SQF states when they get 
closer. For this purpose, the dot is required to have suffi- 
ciently high magnetization. We then expect the GVs can 
be stable if the dot's size is small enough and its magneti- 
zation is high enough. These physical arguments are our 
starting point. In this section, we consider a small circu- 
lar dot R < A, placed upon a SC thin film. Previously, 
we studied the vortex states in the absence of antivor- 
tices and found that at first the L = 1 state appears, 
and, with further increase of magnetization, L = 2,L = 3 
states occur in turn. In addition, the number of states 
depends on the magnetization and the dot's size. We 
anticipate the same sitution when the dot's size is small 
enough to make the GV states stable. Contrary to our 
previous results, we expect L = 1, L g = 2, L g = 3 and so 
on, when only GV states are present. The first question 
we ask here is how small does the dot have to be so that 
the L g = n state will be preferred over the L = n state. 
The necessary conditions for this case are A£ < A„ and 
A^ < 0. We first evaluate this problem analytically, and, 



in doing so, we disregard the antivortices to simplify the 
calculations. They will be discussed at the end of this 
section. At the R < A limit, the effective energy of a GV 



state L g = n reads 
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where Ti is the position of ith antivortex from the dot's 
center. This problem can be simplified by using circular 
symmetry in the system. If the vortices are considered to 
be situated on a ring of radius r v > A such that nearest 
vortices are equally distant from each other, the distance 
between ith and jth antivortices can be expressed as |r^ — 
rj| = 2r v \ sin(7r(z — j)/n)\. Thereafter, taking the radius 
of the ring as r v , the sum in the second term of Eq.(8) 
can be replaced by J2k=i( n ~ fc) ln(A/(2r„ sin(7rfc/n))). 
Under this assumption, the last term of Eq.(8) becomes 
(nm$ -R/2A)(l - r 2 /(4i? 2 )). Note that this assumption 
is valid for up to a few vortices. We expect that new rings 
occur as new vortices come out. In the next section, we 
show that the first ring carries 6 vortices on it. However 
this suffices for our analysis, since it is very unlikely to 
see a GV state with vorticity n = 6. As a result, the 
effective energy is now a function of r„. Optimizing the 
effective energy with respect to r v , we find its equilibrium 
value as 



r v = 4/4(71-1) 



eoRA 
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Directly substituting Eq.(9) into Eq.(8), we find 
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where C n = IT^ 1 sin(7r/c/n)"- fc and e = 2.71828. Let 
us compare Eq.(10) with the energy of a system in the 
presence of a GV with vorticity n 
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The above equation gives the upper limit for the dot's 
size as a function of magnetization of vorticity. When 
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the dot's size is smaller than the upper limit, at first the 
L = 1 state appears. The L g = 2 state follows next. 
With a further increase of magnetization, new GV states 
appear spontaneously. This picture is very similar to our 
previous results, except we have the GV states instead of 
the SQF states. In order for the L g — n state to appear, 



the necessary conditions are < A ! ^ l _ 1 and < 0. 
These conditions give us the lower limit as 

J > U • (13) 



(a) L g = 2 vortex state 



a 3 - 



(b) L g = 3 vortex state 



FIG. 2. The shaded regions show where L g = 2 and L g = 3 states can appear. 



As seen from Fig. 2, as the vorticity of the GV states 
increases, the region in which they might appear shrinks. 
This suggests that the high vorticity GV states are more 
unlikely in usual circumstances. As a result of analytical 
calculations, we expect that the giant vortices can occur 
when R/£ ~ 3,4. To get the better result, we do nu- 
merical calculations for small values of R/ A as analytical 
calculations suggest. In our numerical calculations, we 
set k — 15 and d/A = 0.05. The results of the numerical 
calculations are depicted in Fig. 3. We find that the gi- 
ant vortices dominate the vortex states when R/A = 0.3 
= kR/A = 4.5). However, along with the further 
increase of the dot's size to 0.35 ( R/£ = 5.25), we be- 
gin to see L = 2 and L = 3 states. However, we do not 
see L = 4 or other states with a single quantum flux. 



The physical reason for this is that the dot is still small 
enough so that the vortices can find the enough space 
to decay into SQF states. When the size is increased to 
0.4 (R/£ — 6) , we do not see any GV states anymore. 
The phase diagram of GV states is given in Fig. 4. The 
very bottom curve separates the no vortex region ( area 
below the curve ) from the region where the L = 1 state 
becomes stable ( area between the 1st and 2nd curves). 
This curve is calculated from the equation Ai = 0. Next, 
the L g = 2 state follows. The corresponding curve, which 
separates the L = 1 and L g — 2 states ( area between 2nd 
and 3rd curves ), is found from Aj = Ai. The following 
transitions are found from the equation A^ = for 
ra = 2,3.... 



4 




m<J» /e 



(c) R/A = 0.4 




RJX 

FIG. 4. Phase diagram of GV states. 
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Next, we consider the possibility of AV states together 
with GV states. The total flux due to the magnetic dot is 
zero throughout the entire SC film. The flux due to the 
dot causes the spontaneous appearance of GV or SQF 
states with the dot boundaries. By the same token, the 
flux due to the dot outside the dot boundaries might 
cause spontaneous AV states, when the dot's magnetic 
field penetrates into the the film. Since the total flux is 
zero over the entire SC film, we expect then the number 
of antivortices with single quantum flux equals the GV 
states's vorticity. In principle, the total GV's vorticity 
might differ from the number of SQF AV states in such a 
way that the net vorticity is positive. For the sake of sim- 
plicity, we do not consider this case. However, if stable 
AV vortices do not exist, this scenerio is quite unlikely. 
Because of circular symmetry, we assume that they are 
located on a ring of radius R < r a < A. Again, we stress 
that this assumption is valid only for a few vortices. It 
is, however, sufficient to get a qualitative idea about AV 
states. Under the assumption that AV states arc located 
within the region of A (r a < A limit), the effective energy 
of a GV state L g — n with n antivortices surrounding it 
reads 

A r ~ (n 2 + n)e v + n(n - l)e In — - 2n 2 e a In — 

2r„ r„. 



nm$ R /, R 
-2e lxiC„ ^^—fi- — 
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Minimizing the above equation with respect to the ring's 
radius r a , we find 
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Directly substituting Eq.(15) into Eq.(14), we obtain 
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The GV state L g — n is more energetically favorable than 
the L ga = n state, when A£ 9 < A£ and A° 9 < 0. Using 
the initial approximation R < r a < A and Eq.(15), we 
determine the boundaries of to$o/£o value as 

, A to$o , , A 2 , , 

4(n+l)-<^<4(n+l) — . (17) 

Next, we compare A^ a and A^ for several values of vor- 
ticity n within the limits of m<j> /eo which are given in 
Eq.(17). It turns out than the giant vortices are always 
more favorable than the case in which they appear to- 
gether with antivortices in R < r a < A limit. At the 
opposite limit, r a > A , it is quite unlikely to have 
stable states together with antivortices. Since the effec- 
tive energy goes to n(n + l)eo In k — nm$ R/ (2A) when 



r a > > A. The energy difference neo In k can be over- 
come only when R ~ £. However, we do not consider 
that case in this work because of the inaccuracy of Lon- 
don approximation at the small k limit. Actually, when 
Eq.(2) is optimized together with antivortices, it mani- 
fests instability for vortex states, that is, high vorticity 
states become more favorable, and no low vorticity states 
appear. The reason for this is logarithmic divergence in 
the Hq(x) — Yq(x) function. As x — > 0, this function 
goes logarithmically to infinity, which suggests that in 
the London approximation, a vortex-antivortex interac- 
tion is not taken care of accurately when the distance 
between the vortex-antivortex pairs is on the order of a 
few £. For this situation, we believe that the best ap- 
proach is the nonlinear Landau Ginzburg equation. 



III. SQF VORTEX STATES 

In the second section, we found that, when the dot's 
size is as large as about 0.4A (6£), only SQF vortex states 
occur. We devote this section to the analysis of geo- 
metrical patterns formed by SQF states and their phase 
diagram. In this part, we focus on two possible cases; 
vortices with and without antivortices. The former case 
has been previously analyzed for up to three vortices [19]. 
We determined the geometrical configurations of vortices 
in the ground state. Due to symmetry, the centers of the 
two vortices are located on a straight line connecting the 
vortices with the center of the dot at equal distances from 
the center. The occurrence of two vortices can be exper- 
imentally detected as a violation of circular symmetry of 
the field. For three vortices, the equilibrium configura- 
tion is a regular triangle. In this section, we consider fur- 
ther cases. Due to the circular symmetry, vortices seem 
to favor location on a circle such that they are equally dis- 
tant from nearest neighbors. However, there is another 
possible case in which a vortex is situated at the dot's 
center while others are located on a circle as described 
just above, in the absence of antivortices. 

In order for N vortices to appear, the necessary condi- 
tion is that Ajy < and A at < Ajv-i- Using this criteria, 
we can determine in what configurations and order the 
vortices appear. To this end, we study only vortices with 
positive vorticity that are situated under the dot. The 
next step is to minimize with respect to the positions of 
vortices. We first start with one vortex. It turns out 
that it appears at the center of the dot. Ai is a func- 
tion of two dimensionless parameters mi'o/eo and R/A. 
Ai defines a critical curve that separates regions with or 
without vortices, as is depicted in Fig. 5. Next, we cal- 
culate A2 for two vortices. Our calculations show that 
they are located on a straight line connecting the vor- 
tices with the center of the dot at equal distances from 
the center. The equation A2 = Ai gives the second curve 
in Fig. 5, which separates the regions where one vortex 
and two vortices appear. That is, the area between the 
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first curve and second curve indicates where one vortex 
occurs, while the region between the second curve and 
third curve is for two vortices. The equilibrium config- 
uration of three vortices is a regular triangle. A further 
increase of m<5>o/eo makes other vortex states more en- 
ergetically favorable. Up to seven vortices tend to locate 
on a ring and equally distant from each other. Namely, 
four vortices form a square, whereas five vortices appear 
at the corners of a pentagon. The equilibrium configura- 
tion for six vortices is a hexagon. However, the seventh 
vortex appears at the dot's center, while the other six 
vortices form a hexagon around it. Geometrical patterns 
of up to seven vortices are depicted in Fig. 6. In prin- 
ciple, there exists an infinite series of such transitions. 
Here, we limit ourselves to seven vortices. The curves in 
a phase diagram of these seven vortices are obtained from 
An = Ajv-i equation, which is a function of R/A and 
m$o/^o. Fixing the former, we calculate the latter vari- 
able. We do this for various values of R/A. Fitting the 
points that are obtained from numerical calculations, we 
find the generic function m$ /eo w 4.22 N/ [R/K) u where 



v = 0.99 ±0.03. 




FIG. 5. Phase diagram of the vortex states. 
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(a) One vortex 



(b) Two vortices 



(c) Three vortices 



(d) Four vortices 
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(e) Five vortices 



(f) Six vortices (g) Seven vortices 



FIG. 6. Vortex states 



Equilibrium positions of SQF states change signifi- 
cantly with increasing m$o/eo- That is, the larger 
to$o/£o: the more vortices are pushed towards the dot's 
center. 

Next, we study SQF states in the presence of AV states. 



First, we optimize Eq.(2) with respect to positions of vor- 
tices and antivortices. In our numerical calculation, we 
start from, R/A = 0.4 and and try larger values of R/A. 
We limit ourselves to the L a = 2 and L a = 3 states. Equi- 
librium configurations of L a = 2 and L a = 3 states are 
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depicted in Fig. 7. These configurations are independent 
of the dot's size, wheras equilibrium positions are depen- 
dent on m$ /£o- The further increase of this parameter 
causes the further push of vortices towards the dot's cen- 
ter, and antivortices away from the dot's boundaries. As 
seen in Fig. 7, antivortices are aligned with vortices and 
appear on a ring outside the dot. Next, we compare the 



effective energies of these states with SQF states with- 
out antivortices around them. The comparison of effec- 
tive energies with and without AV states shows that AV 
states are not energetically favorable. However, as dis- 
cussed in the previous section, for dots whose size is on 
the order of a few £, there might be stable AV states. 
However, the London theory does fail in that case. 





(a) L a = 2 vortex state (b) L a = 3 

FIG. 7. Equilibrium configurations of L a = 2 and L a = 3 states. 



IV. CONCLUSIONS 



We studied the vortex states that occur spontaneously 
in the ground states of a SC thin film with a single FM 
dot grown upon it, in the London approximation. Our 
calculations showed that GV states are more energeti- 
cally favorable when the dot's size is up to 0.3A ( 4.5£ ) 
. Between R = 0.3A ( 4.5£ ) and R = 0.4A ( 6£ ) , both 
GV and SQF states appear. Beyond 0.4A ( 6£ ), only 
SQF states are stable. 

Furthermore, we determine the geometric patterns of 
up to seven vortices. Our results show that vortices form 
regular geometric patterns, such as a regular triangle, 
square, pentagon and hexagon. However, in the case 
of seven vortices, one vortex occurs at the dots center, 
while other six sit at the corners of a hexagon. The 
phase diagram of SQF vortex states obeys an equation 
N/(R/A), which suggests that high-vorticity states pre- 
fer larger dots. 

We also studied the cases together with AV states. 
However, our calculations did not show any stable vortex 
states together with antivortices. This result might seem 
implausable since one would expect the presence of AV 
states due to the zero flux condition over the entire film. 
However, the zero flux condition can also be satisfied by 
the presence of a supercurrent outside the dot bound- 
aries, which circulates in the opposite direction to the 
currents created by vortices under the dot. Furthermore, 
the London approach does not enable us to analyze the 
vortex and antivortex states when the dot's size is as large 
as a few £, because it fails for distances on the order of 
coherence length. This failure is significantly transparant 
in the presence of AV states. However, positive vorticity 
states can be studied through the London approach. 

In closing, we studied SQF vortex , GV and AV states 



due to a magnetic dot in the London approximation. Ex- 
perimentally, these states can be checked by magnetic 
force microscopy techniques. 
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